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The cyclic variation of solar activity is both irregular and intermittent. We have sought to 
isolate and illuminate the physical mechanisms of this behavior and to provide a mathematical 
description of it. Our work has brought out three ingredients of the solar cycle that we believe 
to be central to its operation. (1) The seat of the solar cycle is in a shear layer just below the 
solar convection zone. We have investigated the structure of this layer (which we call the solar 
tachocline) in some detail. (2) The spatio-temporal development of the solar cycle is represented 
by the propagation of robust solitary waves which are affected by dissipation and instability. We 
have studied the structure and interactions of such waves, which we call solitoids. (3) On top of the 
simple propagative behavior of the solar solitoids there are intermissions during which the number 
of sunspots remains quite small. We attribute these intermissions (such as the Maunder minimum) 
to a form of interaction between the convection zone and the tachocline which is characteristic of 
a process that we have developed and that we call on/off intermittency. These three ingredients 
make up some of the key features of the solar cycle and may be expected to play a role in future 
simulations of the solar cycle, 
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CHAOTIC DYNAMICS OF THE SOLAR CYCLE 

February 1996 


The Goal 


The existence of sunspots has been known for millenia. It has been known for 
sixty years that spots on the sun occur when strong magnetic flux tubes protrude 
from the solar surface. The strong fields locally inhibit convection and so reduce 
the energy flow from the sun. The region of reduced heat flux naturally is darker 
than its surroundings and is called a sunspot. 

In the middle of the last century it was discovered that there is a cyclic variation 
of the sunspot number. This variation consists of a rise and fall of the number of 
spots on the sun in an 11.1-year cycle which is not strictly periodic. Moreover, 
the location of the principal solar activity, as demarked by sunspots, is localized in 
space and time to a band of latitudes that drifts from midlatitude to the equator 
over the course of the solar cycle. This spatio-temporal variation is expressed in the 
so-called butterfly diagram. 

The aim of our project has been to isolate the essential physical causes of the 
cyclic behavior of spot formation and to pursue their mathematical and physical 
ramifications for the sun, for other stars and in the more general context of nonlinear 
science. The main results obtained in this effort are these: (1) A localization of the 
seat of solar activity to a layer just beneath the solar convection zone. We call 
this layer the solar • tachocline. (2) The explanation of the butterfly diagram of 
solar activity in space and time in terms of nonlinear solitary waves (which we call 
solitoids) and the development of the theory of such waves in several directions. 
(3) An explanation of the dearth of sunspots for seventy-five years in the time of 
Louis XIV by a mathematical mechanism we have called on/off intermittency and 
the connection of this process to the solar tachocline by a simple dynamical model. 

In this report, we briefly outline these three developments. We also add a 
bibliography below that refers to the several other projects that grew out of these 
studies. A more extensive summary is in preparation for presentation in the Enrico 
Fermi School of Physics this June and will be forwarded at that time. 



The Solar Tachocline 


The fact that solar activity is, at any time, restricted to a fairly narrow band 
of latitude has led us to develop the idea that the source of sunspots is confined 
to a layer in the sun of comparable, if somewhat smaller, depth than the width 
of the excitability band. The convection zone of the sun occupies the outer third 
of the solar radius and is therefore too deep to meet this requirement. We have 
therefore developed the idea that sunspots form in a layer just below the convection 
zone. The existence of such a layer, adumbrated twenty-five years ago, has now 
been confirmed by helioseismology and we have called it the solar tachocline. It 
mediates the transition from the outer, strong differential rotation to the more rigid 
rotational behavior of the inner sun. 

In our picture, the active solar convective layer, which occupies the outer third 
of the solar radius, is the analogue of the earth’s atmosphere while the inner, deep 
core of the sun corresponds to the earth’s oceans. In this analogy, the tachocline 
corresponds to the oceanic thermocline, and this is the source of its name. We have 
made a preliminary theory of its structure which is in fair agreement with the best 
results of helioseismology. The solar cycle is then a sort of magnetic weather system 
of the sun in which magnetic rainfall descends from the convection zone only to be 
returned from the tachocline in the form of magnetic storms — the sunspots. 



Fig. 1 . Structure of the Solar Tachocline (from “The Solar Tachocline” by E.A. 
Spiegel and J.P. Zahn, Astron. & Astrophys., 265, 106-114, 1993). 
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It now seems to be generally accepted that the failure of early simulations 
directed at reproducing the solar cycle accurately may be due to the omission of a 
tachocline from the modeling. In the wake of recent highly successful simulations 
of the geodynamo by Glatzmaier and Roberts we may hope that the solar problem 
may soon be similarly studied and that the inclusion of a tachocline will be a key 
ingredient just as inclusion of a solid core was essential in the terrestrial work. 


Solar Solitoids 


The plot of the latitude of solar activity as a function of time that is called 
the butterfly diagram shows that solar activity occurs in a limited range of latitude 
that drifts during the solar cycle from midlatitude to the equator. There is similar 
activity in each hemisphere, though in detail there is significant asymmetry about 
the equator. 

Our explanation of the butterfly pattern is that there is a pair of solitary waves, 
one in each hemisphere, each propagating toward the equator. These nonlinear 
waves arise because of an a magnetohydrodynamic instability of the tachocline that 
develops as plunging plumes in the convection zone bring magnetic field down to be 
intensified by the shear of the tachocline. It is in the nature of mhd instabilities that 
they often occur as overstability or growing oscillations. Whatever the details of the 
mechanism, its nonlinear equation is generically described by a nonlinear equation 
known as the complex Ginzburg-Landau equation. We have adapted this equation 
to the case of the tachocline by allowing for latitude dependence of its governing 
parameters. In the examples of this form of equation that we have studied, solitary 
waves of the right kind are produced and propagate in a convincing manner toward 
the equator. 



Fig. 2. Propagation of a Solar Solitoid in the Northern Hemisphere (from “Waves 
of Solar Activity,’ by M.R.E. Proctor and E.A. Spiegel, in The Sun and Cool Stars, 
(Springer-Verlag) 1991). 
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Like solitons that are found in nondissipative integrable systems these solitary 
waves possess a good deal of stability. Unlike solitons they are subject to dissipation 
and instability. These extra features make them more interesting in that their 
parameters (amplitude, breadth, velocity) are not constants, as for solitons, but are 
subject to variations that may be complicated. In the cases where the dissipation 
and instabilities are in a reasonable balance on average, the parameters of these 
solitary waves are governed by dynamical systems which control the behavior of the 
solitary waves. We call these such solitary waves solitoids. 

We have derived equations of motion (dynamical systems) for the solar solitoids 
that propagate in the tachocline. When we allow for the latitudinal variations of 
the properties of the tachocline, we find a reasonable description of the butterfly 
diagram. The situation is enriched when the solitoids in the two solar hemisphere 
are allowed to interact. In this case, we find mildly irregular behavior and north- 
south asymmetry as one sees in the observed solar cycle. 

To clarify these aspects of the problem we have also made detailed studies of 
the interactions of solitary waves in general and of the spatio-temporal chaos this 
produces. The latter theory has a variety of applications in applied mathematics 
and we have developed the theory of interacting solitoids for such studies, with 
applications to biomathematics and other branches of astrophysics. 


On/Off Intermittency 


For nearly three hundred years, the number of spots on the sun at any time 
has fluctuated up and down on the average time of 11.1 years. This variation 
resembles that of a chaotic oscillator such as those third-order systems that have 
been much studied in recent decades. But, during the time of Newton, sunspots 
became quite rare and this lowering of activity has puzzled investigators ever since. 
It is especially intriguing because there are hints that this dearth of sunspots, known 
as the Maunder minimum, is correlated with little ice ages. 

We have devised a mechanism to explain the near disappearance of sunspots, 
one that fits in with our view that the solar cycle originates in the tachocline. The 
mechanism, called on/off intermittency, is a general one that explains how systems 
alternate between high activity of a certain kind and quiet periods of inactivity. 
The alternation of a turbulent fluid between turbulent and laminar phases is an 
example of this kind of behavior. 

A system that exhibits on/off intermittency has two components, an observable 
oscillator that is activated by a driver that is a larger system whose variables are 
neither predictable nor observable. The sudden outbursts of activity of the oscillator 
cannot be predicted either. In the application of this mechanism to the solar cycle 
the tachocline is the observable oscillator and the solar convection zone is the driver. 
The magnetic field that rains down from the convection zone enters the tachocline 
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where it is twisted into a strong toroidal field by the characteristically strong shear 
across that layer. Since the density is stably stratified in the tachocline, the field 
is confined there until it becomes strong enough to erupt into the convection zone 
where it is buoyed up to form sunspots on the surface. If the level of magnetic rainfall 
drops below a certain critical level, the instability of the tachocline does not occur 
and the main solar cycle ceases, though dynamo activity in the solar convection 
zone continues. This picture gives rise to the process of on/off intermittency in the 
sun and is responsible for intermissions of the cycle. 
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Publications 


Each of these ingredients of the solar cycle has been the object of separate 
lines of investigation and it will need a few more years to work out the details of 
their intermeshing. We give now a list of papers in which these ideas and related 
developments appear: 


“Waves of solar Activity,” (M.R.E. Proctor and E.A. Spiegel), in the sun and the 
Cool Stars, D. Moss, G. Ruddinger and I. Tuominen, eds. (Springer-Verlag) 
117-128,1991. 

“Interacting Localized Structures with Galilean Invariance,” (C. Elphick, O. Regev, 
G.R. Ierley and E.A. Spiegel), Phys. Rev. A, 44, 1110-1122, 1991. 

“On Thermonuclear Convection: I Stability Theory,” (S. Ghosal and E.A. Spiegel), 
Geophys. & Astrophys. Fluid Dyn., 61, 161-179, 1991. 

“Vortices on accretion disks,” (with M.A. Abramowicz, A. Lanza, E.A. Spiegel and 
E. Szuszkiewicz), Nature, 356, 41, 1991. 

“Non-linear instability of viscous plane Couette flow, I. Analytical approach to a 
necessary condition,” (B. Dubrulle and J.-P. Zahn) J. Fluid Mech., 231, 561, 
1991. 

“Convective penetration in stellar interiors,” (B. Dubrulle and J.-P Zahn), Astron. 
Astrophys., 252, 179, 1991. 

“Convective penetration into stellar radiation zones,” (J.-P. Zahn) in Challenges 
to Theories of the Structure of Moderate Mass Stars, D.O. Gough and J. 
Toomre, eds. (Springer) 225, 1991. 

“Effects of turbulence on large scale motions: towards a model of coherent struc¬ 
tures,” (M. Rieutord and J.-P Zahn), in The Sun and Cool Stars: Activity, 
Magnetism and Dynamo. I. Tuominen, ed. (Springer) 33, 1991. 

“Evidence for convective penetration in the Sun,” (G. Berthomieu, P. Morel, J. 
Provost and J.-P. Zahn), in 7th Cambridge Conf. on Cool Stars, Stellar 
Systems and the Sun, M. Giampapa, ed. (Conf. Proc. Astron. Soc. Pacific) 
240, 1991. 

“Rotating stars,” (J.-P. Zahn) in Lecture notes prepared by students for the Pro- 
ceeedings of the Summer Study Program in G.F.D., R. Salmon, ed. (Woods 
Hole Oceanogr. Inst. Tech. Rep.; WHOI 91-03) 52, 1991. 

“Convection and turbulence in Stars,” (J.-P. Zahn) in Late Stages of Stellar Evo¬ 
lution, C.B. de Loore, ed. (Springer) 355, 1991. 
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“Convection-induced shears with general planforms,” (J. Massaguer, E.A. Spiegel 
and J.-P. Zahn) Phys. Fluids A, 4, 1333, 1992. 

“The turbulent tachyline,” (E.A. Spiegel and J.-P. Zahn), in 7th Cambridge Conf. 
on Cool Stars, Stellar Systems and the Sun, Astron. Soc. Pac. Conf. Series 
26, 158, 1992. 

“The solar tachocline,” (E.A. Spiegel and J.-P. Zahn) Astron. <fc Astrophys., 265, 
106-114, 1992. 

“Effect of horizontal turbulent diffusion on the transport by meridional circulation,” 
(B. Chaboyer and J.-P. Zahn), Astron. Astrophys., 253, 173, 1992. 

“Convection-induced shears of general planforms,” (J. Massaguer, E.A. Spiegel, and 
J.-P. Zahn) Phys. Fluids, A 4 (7), 1333, 1992. 

“Circulation and turbulence in rotating stars,” (J.-P. Zahn), Astron. Astrophys., 
265, 115, 1992. 

“On-off intermittency: A mechanism for bursting” (N. Platt, E.A. Spiegel and C. 
Tresser) Phys. Rev. Lett., 70, 279-282, 1993. 

“Equilibria of rapidly rotating polytropes,” (N.J. Balmforth, E.A. Spiegel and L.N. 
Howard) Mon. Not. Roy. a str. Soc., 260, 253-272, 1993. 

“Astrophysical Fluid Dynamics,” (E.A. Spiegel) in Astrophysical Fluid Dynamics, 
J.-P. Zahn and J. Zinn-Justin, eds., Les Houches Session XLVII, (Elsevier 
Science Pubs.) 3-33, 1993. 

“Patterns of Aperiodic Pulsation,” (E.A. Spiegel) Astrophys. & Space Sci., 210, 
33-49, 1993. 

“The Intermittent Solar Cycle” (N. Platt, E.A. Spiegel and C. Tresser) Geophys. 
& Astrophys. Fluid Dyn., 73, 146-161, 1993. 

“Turbulence in stellar interiors,” (J.-P. Zahn) in Inside the stars, A. Baglin and W. 
Weiss, eds. (Conf. Proc. Publ. Astron. Soc. Pacific) 236, 1993. 

“Instabilities and turbulence in rotating stars,” (J.-P. Zahn), in Les Houches, session 
XLVII : Dynamique des Fluides Astrophysiques, J.-P. Zahn and J. Zinn- 
Justin, eds. (Elsevier Science Publ.) 561, 1993. 

“Seismological constraints on convective penetration in the Sun,” (G. Berthomieu, 
P. Morel, J. Provost, and J.-P. Zahn) in Inside the stars, A. Baglin and W. 
Weiss, eds. (Conf. Proc. Publ. Astron. Soc. Pacific) 60, 1993. 

“Topological entropy of one-dimensional maps: approximations and bounds” (N.J. 
Balmforth, E.A. Spiegel and C. Tresser) Phys. Rev. Lett., 72, 80-83, 1994. 




“Advection of Vector Fields by Chaotic Flows” (N.J. Balmforth, P. Cvitanovic, G.R. 
Ierley, E.A. Spiegel and G. Vatay) in Stochastic Processes in Astrophysics, 
J.R. Buchler and H.E. Kandrup, eds. Ann. N.Y. Acad. Sci., 706 148-160, 
1994. 

“Autogravity Waves in a Polytropic Layer” (Z.S. Qian and E.A. Spiegel) Geophys. 
& Astrophys. Fluid Dyn., 74, 225-244, 1994. 

“Instability in extended systems,” (L. Tao and E.A. Spiegel), in Proceedings of the 
Summer Study Program in G.F.D., Woods Hole Oceanographic Institution, 
R. Salmon, dir., pp. 177-184, WHOI-94-12, 1994. 

“Chaotic Pulse Trains,” (N.J. Balmforth and G.R. Ierley and E.A. Spiegel) SIAM 
J. Applied Math., 54, 1291-1234, 1994. 

“The Chaotic Solar Cycle,” (E.A. Spiegel) in Lectures on Solar and Planetary Dy¬ 
namos, M.R.E. Proctor and A.D. Gilbert, eds., Pubs. Newton Inst. (Cam¬ 
bridge Univ. Press), 245-266, 1994. 

“Penetration below a convection zone,” (N.E. Hurlburt, J. Toomre, J. Massaguer, 
J.-P. Zahn), Astrophys. J., 421, 245, 1994. 

“Mixing processes and stellar evolution,” (J.-P. Zahn), Space Science Reviews, 66, 
285, 1994. 

“Plumes in stellar convective envelopes,” (M. Rieutord and J.-P. Zahn), Astron. 
Astrophys., 296, 127, 1995. 

“Instability of Rapidly Rotating Polytropes,” (N.J. Balmforth, L.N. Howard and 
E.A. Spiegel) SIAM J. Applied Math., 55, 298-331, 1995. 

“Checkerboard Maps,” (N.J. Balmforth, E.A. Spiegel and C. Tresser) Chaos, 5, 
216-226, 1995. 

“A Prelude to Stellar Convection Theory,” (E.A. Spiegel) in Physical Processes 
in Astrophysics, I.W. Roxburgh and J.-L. Masnou, eds., Lecture Notes in 
Physics ( Springer-Verlag) 129-152, 1995. 
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